home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / DSPDI.z / DSPDI
Text File  |  1998-10-30  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDSSSSPPPPDDDDIIII((((3333FFFF))))                                                            DDDDSSSSPPPPDDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DSPDI   - DSPDI computes the determinant, inertia and inverse of a double
  10.      precision symmetric matrix using the factors from DSPFA, where the matrix
  11.      is stored in packed form.
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.       SUBROUTINE DSPDI(AP,N,KPVT,DET,INERT,WORK,JOB)
  16.  
  17. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  18.      On Entry
  19.  
  20.      AAAAPPPP DOUBLE PRECISION (N*(N+1)/2)
  21.         the output from DSPFA.
  22.  
  23.      NNNN INTEGER
  24.         the order of the matrix A.
  25.  
  26.      KKKKPPPPVVVVTTTT INTEGER(N)
  27.         the pivot vector from DSPFA.
  28.  
  29.      WWWWOOOORRRRKKKK DOUBLE PRECISION(N)
  30.         work vector.  Contents ignored.
  31.  
  32.      JJJJOOOOBBBB INTEGER
  33.         JOB has the decimal expansion  ABC  where
  34.         if  C .NE. 0, the inverse is computed,
  35.         if  B .NE. 0, the determinant is computed,
  36.         if  A .NE. 0, the inertia is computed.
  37.         For example, JOB = 111  gives all three.  On Return Variables not
  38.      requested by JOB are not used.
  39.  
  40.      AAAAPPPP contains the upper triangle of the inverse of
  41.         the original matrix, stored in packed form.
  42.         The columns of the upper triangle are stored
  43.         sequentially in a one-dimensional array.
  44.  
  45.      DDDDEEEETTTT DOUBLE PRECISION(2)
  46.         determinant of original matrix.
  47.         DETERMINANT = DET(1) * 10.0**DET(2)
  48.         with 1.0 .LE. DABS(DET(1)) .LT. 10.0
  49.         or DET(1) = 0.0.
  50.  
  51.      IIIINNNNEEEERRRRTTTT INTEGER(3)
  52.         the inertia of the original matrix.
  53.         INERT(1)  =  number of positive eigenvalues.
  54.         INERT(2)  =  number of negative eigenvalues.
  55.         INERT(3)  =  number of zero eigenvalues.  Error Condition
  56.  
  57.      AAAA division by zero will occur if the inverse is requested and  DSPCO  has
  58.      set RCOND .EQ. 0.0 or  DSPFA  has set  INFO .NE. 0 .  LINPACK.  This
  59.      version dated 08/14/78 .  James Bunch, Univ. Calif. San Diego, Argonne
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDSSSSPPPPDDDDIIII((((3333FFFF))))                                                            DDDDSSSSPPPPDDDDIIII((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      Nat. Lab.  Subroutines and Functions BLAS DAXPY,DCOPY,DDOT,DSWAP Fortran
  75.      DABS,IABS,MOD
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.